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4.6 Selective withdrawl of thermally stratified fluid 
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We now extend the analysis in the last section and consider the slow and steady flow of 
a thermally stratified fluid into a two- dimensional line sink. 




Pigur* 4.6. i : Sketch of velocity prafSlos st-msa th^ fa/er dralrlng Into a Sine sSnkr 
(AdBpted tram Katt. J. Ftuld Macnanlca, 24. 



Thermal diffusion and convection now comes into play. 
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4.6.1 Governing equations 

We begin with the general law of mass conservation 

^ + V ■ (pg) = ^ + q- Vp + pV- (4.6.1) 

In environmental problems the range of temperature variation is within a few tens of degrees. 
The fluid density varies very little and obeys the following equation of state 

p = Po[l-(3iT-T,)] (4.6.2) 

where T denotes the temperature and /? the coefficient of thermal expansion which is usually 
very small. Hence 



pV-^* \ p 
and 

^ ^ ^ « 1 
Vu p 

It follows that (4.6.1) is well approximated by 

V-g = (4.6.3) 

which means that water is essentially incompressible. In two dimensions, we have 

«x + = (4.6.4) 

Next, energy conservation requires that 

dT 

— + q-VT^DV^T (4.6.5) 
ot 

Let 

T = T + T' (4.6.6) 

where T represents the static temperature when there is no motion, and T' the motion- 
induced temperature variation. Therefore, 

T - r„ = (T(z) - To) + T'{x, z, t) (4.6.7) 

and 

&r' - — 

— + q-VT + q-VT'^DV^T + DV^T' (4.6.8) 
ot 

The static temperature must satisfy 

V^T = (4.6.9) 
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In a large lake with depth much smaller than the horizontal extent, the static temperature 
is essentially uniform horizontally. The Laplace equation reduces to 

(PT dT 
D—— = 0, implying — — = constant (4.6.10) 
az^ dz 

The dynamic part is then gorvened by 

dT' dT' dT' dT 

The exact equations for momentum balance are, in two dimensions, 

pff + g-V^)=-|+/.V^« (4.6.12) 



where p denotes the static part, which must satisfy 

= - gpo [l - /3(r - TJ] (4.6.14) 

Taking the differenece of the two preceding equations, we find the equation for the dynamic 
part 

4.6.2 Approximation for slow and steady flow 

For sufficiently slow flows, inertia terms can be ignored. Expecting that vertical motion 
is suppressed, we further assume that the vertical length scale 5 is much smaller than the 
horizontal scale L, so that d/dx <C d/dz. The 2-D momentum equations can then be 
simplified to 

= -^^ + gPpX + (4.6.17) 



Similarly we can linearize (4.6.11) to get 



dT d'^T' 

-T.-'^l^ (4.6.18) 



Together (4.6.4), (4.6.18), (4.6.16) and (4.6.17) complete the lineaized governing equations. 
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Eliminating p from (4.6.16) and (4.6.17), we get 



Since 



d\ or 



- = O - « 1, ^ = O - I « 1 
u \L u. \L, 



(4.6.19) 



we can omit the second term on the left of (4.6.19). In terms of the stream function defined 

by 

n = V., w = -iJx (4.6.20) 

(4.6.19) becomes 



Equation (4.6.18) can be written as 





gPPodT' 




jjL dx 


as 


dT 




az 


dz^ 



(4.6.21) 



(4.6.22) 



We now have just two equations for two unknowns and T'. The boundary conditions are 

T'u, u; i 0, as z\ ±oo (4.6.23) 



or 



V'.V^TjO, as 0T±oo. (4.6.24) 
Let the volume rate of withdrawal be prescribed, we must then require the integal condition: 



u dz — —q, implying ijj{x, z — oo) — ip{x, z — — oo) = q. 



(4.6.25) 



4.6.3 Normalization 

Let 



^^q^*, T^ToT*, x^Lx*, z ^ 5z* 



(4.6.26) 



Physically it is natural to choose the characteristic depth of thermal gradient as the global 
length scale L: 



L^- 13 



dz 



(4.6.27) 
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^*(oo) - V*(-oo) = 1 



(4.6.34) 
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4.6.4 Similarity solution 

Let us try a one-parameter similarity transformation 



x = X''x, z = X^z, ^ = \''ip, T' = X'^f 



(4.6.35) 



The exponents a, b, c and d will be chosen so that the boundary value problem is formally 
the same as the original one To achieve invariance of (4.6.34), we set c = 0. In addition we 
set 

for (4.6.29), and 

for (4.6.31). Hence, d — a — —46 and a — 2b — —d implying 

b = -d, a = 3b = -3d. (4.6.36) 
These relationships among the exponents suggest the following new similarity variables: 

. HO 



with 



i> = f{0, 



.1/3 



r-1/3 



(4.6.37) 



(4.6.38) 



It is easily verified that these variables are invareiant under the similarity transformation. 
Carrying out the differentiations 



T ^ hi 



a.l/3 ' ^4/3 

1 / 1 ^ \ ^/ 1 



3x4/3 



3; x4/3 



-rh'— - — 
3^ x4/3 + 3x4/3 



we get from (4.6.29) 



(4.6.39) 



Since 



h' 



3.2/3 

z 



X 

1 
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we get from (4.6.31) 



X 2) X 



or 



h" - = 



(4.6.40) 



The boundary conditions are transformed to 



/(oo)-/(-oo) = -l 



and 



as 



C^±oo 



(4.6.41) 



(4.6.42) 



Mathematically, the similarity transformation has enabled us to reduce the boundary 
value problem involving partial differential equations to one with ordinary differential equa- 
tions (4.6.39), (4.6.40), (4.6.41), and (4.6.42). As long as x and z lie on the parabola 

z = const x^/'^, t/j* and T*x*^^^ are the same. From the transformation, we can also deduce 
that the boundary of the zone affected by the flow is a parabola. 



Along the centerline z — ( = 0, the velocity varies as 



Umax ~ V'x ~ a;" 



-1/3 



and the temperature varies as 



^max 



4/3 



(4.6.43) 



(4.6.44) 



(4.6.45) 



The boundary value problem can now be solved by numerical means (such as Runge- 
Kutta). Numerical results by Koh (1966, Fig. 4) are shown in Figure (4.6.2). 

In Koh (1966), stratification in fluid density is associated with the variation of concen- 
tration of a diffusive substance instead of temperature. The fluid density is governed by 
a diffusion equation formally the same as that for temperature here. To use his numerical 
results, /o,/io in his plots are replaced by our f,—h shown here. Extensive discussion on 
experimental confirmation as well as the three dimensional theory for a point sink can be 
found in Koh. 



